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Abstract 



We calculate the spin-drag transresistivity p^ (T) in a two-dimensional elec- 
tron gas at temperature T in the random phase approximation. In the low- 
temperature regime we show that, at variance with the three-dimensional 
low-temperature result [p^[(T) ~ T 2 ], the spin transresistivity of a two- 
dimensional spin unpolarized electron gas has the form p-m(T) ~ T 2 lnT. 
In the spin-polarized case the familiar form p^(T) = AT 2 is recovered, 
but the constant of proportionality A diverges logarithmically as the spin- 
polarization tends to zero. In the high-temperature regime we obtain p+ 1 (T) = 
— (h/e 2 )(TT 2 Ry*/kBT) (where Ry* is the effective Rydberg energy) indepen- 
dent of the density. Again, this differs from the three-dimensional result, 
which has a logarithmic dependence on the density. Two important differ- 
ences between the spin-drag transresistivity and the ordinary Coulomb drag 
transresistivity are pointed out: (i) The InT singularity at low temperature is 
smaller, in the Coulomb drag case, by a factor e~ 4kpd where kp is the Fermi 
wave vector and d is the separation between the layers, (ii) The collective 
mode contribution to the spin-drag transresistivity is negligible at all temper- 
atures. Moreover the spin drag effect is, for comparable parameters, larger 
than the ordinary Coulomb drag effect. 

72.25.-b, 72.10.-d, 72.25.Dc 
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I. INTRODUCTION 



The problem of transporting an electronic spin polarization from a place to another by 
means of an electrical current has attracted tremendous interest in recent years, both in the- 
oretical and experimental circles [1]. A particularly challenging problem is that of injecting 
spin-polarized electrons from a ferromagnet into an ordinary non-magnetic semiconductor 
such as GaAs. This is by no means easy, due to the large impedance mismatch between the 
ferromagnetic injector and the semiconductor [2]. However, the discovery that GaAs, when 
doped with magnetic Mn impurities, becomes ferromagnetic with a Curie temperature T c 
as large as 110K [3] has raised hopes of realizing useful all-semiconductor spin-electronic 
devices in a near future. 

In dealing with spin-polarized currents one must be prepared to treat a situation in 
which electrons of opposite spin orientations travel, on the average, with different speeds. 
For example, the up-spin electrons injected from a ferromagnet into a semiconductor might 
be drifting in one direction, while the down spin electrons remain on the average stationary. 
Whenever such a situation occurs, the existence of Coulomb correlation between up- and 
down-spins becomes important, because it is a source of friction between the two compo- 
nents. In our previous works [4-6] we examined the effect of up-down spin correlations on 
spin polarized transport in a three-dimensional electron gas, and we introduced the concept 
of spin Coulomb drag (SCD) [4]: due to Coulomb scattering, momentum is transferred be- 
tween the spin populations, tending to equalize their average drift velocities. This process 
represents an actual drag exerted by the slower population on the faster and, in the absence 
of an external driving field, leads to a decay of the spin current even in the absence of extrin- 
sic impurities [7] . Later, it was shown that the SCD also limits the diffusion of a spin-packet 
in a doped semiconductor ( [5]) adding to the conventional effect of impurity scattering. 

The quantitative measure of the SCD is the spin-drag transresistivity p-^(T), which is 
defined as the ratio of the gradient of electrochemical potential for up-spin electron (E^) to 
the current of down-spin electrons (ji) when the current of up-spin electrons is zero: 

Ei=p n (T)j h (j T = 0). (1) 

In Ref. [4] we presented a calculation of Pu(T) in a pure three dimensional electron gas 
in the Random Phase Approximation (RPA). The calculation is analogous to that of the 
more familiar Coulomb drag effect [8], but, as we will show momentarily, some important 
differences exist between the two effects. 

The ordinary Coulomb drag (CD) effect involves two distinct quasi-two-dimensional elec- 
tron layers that are spatially separated by an insulating barrier. By contrast, in the SCD, 
we are dealing with a single layer in which the two carrier populations (distinguishable by 
the different spin) share the same space and interact with the same impurities [9]. It follows 
that the SCD can be either a two-dimensional (2D) or a fully three-dimensional (3D) effect, 
depending on the dimension of the space in which the electrons move. We have already 
shown that the 3D spin transresistivity should be measurable in metals [4] and becomes 
larger in semiconductors [5,6], where it can be of the same order of magnitude as the usual 
Drude resistivity. In this paper we concentrate on the 2D case. Due to the absence of spa- 
tial separation between the two spin populations, the bare interaction between parallel and 
antiparallel spin electrons coincides with the interaction between parallel spin electrons, and 
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its Fourier transform is simply v(q) = 2ne 2 /qe where q is the wave vector in the plane and 
e is the appropriate low-frequency dielectric constant [10]. By contrast, in the standard CD 
case, the bare interaction between electrons in the two layers has the 2-dimensional Fourier 
transform v 12 (q) = 2ne 2 e~ qd / qe where d is the distance between the layers. This "innocent" 
difference has two important consequences. First, there is no acoustic-plasmon- mediated 
interaction in SCD, whereas the effect of acoustic plasmons on the temperature dependence 
of the CD in separate layers is quite marked [11]. Second, the effective interaction between 
up and down spins in SCD remains sizable up to values of q comparable to the Fermi mo- 
mentum: the presence of a significant interaction strength at q = 2k F is responsible for a 
characteristic logarithmic singularity in the low-temperature transresistivity. 

We provide analytical expressions for pjj_ in both the low- and high-temperature limits. 
At low temperature we derive analytically the T 2 In T singularity; at high temperature we 
show that becomes density-independent. Both these limits are significant for realistic 
metallic and semiconductor parameters, respectively. We also draw a comparison between 
the spin Coulomb drag and the ordinary Coulomb drag. It is hoped that the theoretical 
results presented in this paper will stimulate experimental work aimed at a direct observation 
of the SCD and a comparison between SCD and CD [12]. 



II. SPIN TRANSRESISTIVITY IN THE 2-DIMENSIONAL ELECTRON LIQUID 

To leading order in the strength of the Coulomb interaction [13], the spin-drag tran- 
sresistivity is given by Eq. (18) of Ref. [4], which, in the static case and after few simple 
rearrangements, can be written as 

Pn(T) = - r / dq- 



e 2 n 1 n i k B T{2ir) D + 1 Jo * D 
ImQoi(q,u)ImQ ol (q,u) 
o U \€(q,u)\ 2 wah 2 (hu}/2k B T)' 

Here k B is the Boltzmann constant, D is the number of spatial dimensions and Cd = 
1, 27r, 4-7T for D — 1, 2, 3 respectively. The dimensionless function Q 0a (q,u) is defined as 

Q a(q, u) = v(q)xoAq, u), (3) 

where Xoa{q, w) is the non- interacting density- density response function (the "Lindhard func- 
tion") for spin a —\ or \ and v(q) is the Fourier transform of the Coulomb interaction. 
e(q,cu) = 1 — Qoi(q,u) — Qoi(q,u) is the RPA dielectric function of the electron liquid, and 
n-f and are the densities of up- and down-spin electrons. In both two and three dimen- 
sions Imxoa can be expressed analytically. In particular, in two-dimensions, starting from 
the familiar expression 

r d 2 k 

Imxo*(q,u) = -7T J j^y n k {5(hu + E k - E k ^ q ) 

- 8{hw - E k + E k _ q )} , (4) 

with E k = h 2 k 2 /2m* and n k = l/[exp((£' fc — fi)/k B T) + 1] (m* is the effective mass), we 
obtain 
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q z a* v 

- F (e^-MI**)]. (5) 

Here a* = h 2 e/me 2 is the effective Bohr radius, k Fa is the Fermi wave vector of u-spin 
electrons, T a = ksT/Ep^ and p, a = p a /E Fa are the temperature and chemical potential of 
a-spin electrons expressed in units of the a-spin Fermi energy E Fa , and 

qv Fa 2k Fa 

where v Fa = Hk Fa / m* is the a-spin Fermi velocity. 
The function F in Eq. 5 is defined as 

, r°° dx 

F z 



r (7) 

Jo ze* 2 + l [) 



In Fig. 1 we present the behavior of Pu(T) calculated from Eq. (2) and D = 2 for three 
different values of the carrier density with the values of the effective mass and dielectric 
constant appropriate for InAs (m* = 0.026m e , m e the free electron mass and e = 13.6). 
As pointed out in [4,6], the spin-transresistivity increases with decreasing densities. As 
a function of temperature, vanishes at T = 0, reaches a maximum about the Fermi 
temperature T F) and decreases for large T. As shown in Fig. 1, the overall scale of the 
effect, at its largest, is set by h/e 2 ~ A.lkQ. This is quite a large value and should definitely 
be observable. 

We now derive and discuss the degenerate and non degenerate limits of Eq. (2). 



III. THE DEGENERATE LIMIT 



In the limit T <^ T F the calculation of Pn(T) can be carried out analytically. The 
temperature dependence of p^, in this regime, is entirely determined by the denominator 
sinh 2 (hco/ '2k bT), which restricts the integral in Eq. (2) to frequencies of the order of ksT. 
We can therefore neglect the temperature dependence of ImQ 0(T and replace this function 
by its zero-temperature limit 



ImQ a{q,uj) 



k 



Fa 



q 2 a* 



e(i-^)\A-e 



-e(i-v 2 )Ji 



(8) 



The behavior of ImQ 0a (q,uj) for small u depends crucially on the value of q. As shown in 
the upper panel of Fig. 2, there are three distinct regions in which ImQo a (q,u) is different 
from zero: 



( J ) *** = k ^ [- 1 + V^ + ^d < q < k Fa [l - v^=g= 

k Fa 



= qBa- 



ImQ 0a (q,uj) 



q 2 a* V k Fa v Fa ' 



(9) 
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(II) q Ba = k Fa 1 - yjl 



2u) 



ImQ 0(7 (q,uj) 



< q < k Fa 



1 + a/1 



2u) 



kFtrVptj 

k Fa uj/(k Fa v Fa ) 
Jl-(q/2k Fa f 



qca 



(10) 



(III) q c , = k Fa [l + v/l-^j < q < k Fa [l + V / TT^g= 

^(7 



q 2 a* V k Fa v Fa 



+ 1 



(11) 



Let us first consider the paramagnetic case = k F ^ = k F ± etc.. In region I the 
product IrnQo^IrnQoi is proportional to u, while the integration region in q space is a shell 
of thickness proportional to uo 2 . The net q-dependence of the integrand (taking into account 
the fact that e(q,0) ~ 1 + 2/(qa*), |e(g,0)| 2 ~ l/(qa*) 2 ) is oc q, giving an extra factor u/v F . 
Since %oo ~ k B T, the contribution of this region to is at least of order T 4 and will be 
disregarded hereafter. 

In region II the product ImQo^ImQoi is proportional to u 2 , while the integration region 
in q space is a shell of thickness ~ 2k F . This would give a contribution of order T 2 (as in 3 
dimensions) were it not for the square root divergence of ImQ Q(T when q approaches 2k F [see 
Eq. (10)]. Due to the "piling-up" of two square-root singularities, one gets an additional 
ln(a>) from the upper limit (q ~ 2k F ) of the integration interval. Thus the total contribution 
of this region to pn is of order T 2 lnT. Notice that this logarithmic contribution would be 
severely suppressed in the case of the ordinary Coulomb drag, since the Fourier transform of 
the Coulomb interaction between electrons in different layers decreases, with increasing q, 
as e~ qd . 

Finally, in region III, the product ImQQ^ImQ^ is again proportional to u>, but the 
integration region in q space is now a shell of thickness ~ ou/v F centered about 2k F . The 
contribution of this region is therefore of order T 2 . 

Combining the contributions of regions II and III together, after some calculations, we 
arrive at the following expression: 



Pn(T) 



t<t f h 2ir r 



2^2 



'k B T 



/( 



e 2 3{V2 + r s ) 2 
k B T 



E, 



T2> 



In 



E f 



+ ln2 + 1 - 



I/ 2 hi ?/ 
it 2 Jo sinh 2 y 



dy 



(12) 



where r s = l/y/nna* is the dimensionless Wigner-Seitz radius , Jq 00 (y 2 In y/ sinh 2 y)dy = 



(13) 



—0.55981 and the function f(x) is defined as 

1 



x 



l + ^±lln(l±^) 



x 



2x 



Thus, in the paramagnetic case the spin trans-resistivity behaves as T 2 InT rather than T 2 . 

This is no longer true if the electron gas is spin-polarized. Let us assume > so 
that fcpf > k F i. For a finite polarization and at sufficiently low temperatures such that 
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k B T <C Ep] and k B T <C E B [, ImQ Qa (q,u) can be approximated with its zero-temperature 
expression and the region of integration in q coincides with the region in which ImQo^q, u) 
differs from zero (see lower panel of Fig. 2). As in the paramagnetic case, it is then possible to 
divide the q- integral into regions in which ImQ 0a takes values given by Eqs. (9)- (11). At the 
upper limit of region II for down spins [Eq. (10)], q = 2/cpj — u/vpi, ImQoi(q,uj) diverges 
as 1/^1 — (q/2k F i) 2 . However, ImQ ^(q,u) remains regular. Therefore, the logarithmic 
divergence, leading to the InT singularity in the transresistivity is cut off. 
Our final result takes the form 

frTiS r«T F h 7r 2 r 2 s (k B T\ 2 
Pn(T) ~ -^^ — ^km -eH 9(r s ,0 (14) 



e 2 3(1 - ef 2 



where the function g is defined as 

9(r s ,0 = I dx- 
Jo 



X 



1 2 



X + 



\J (1 — ax){l — x) 



(15) 



£ = (rij — n^)/(ri| + raj is the degree of spin polarization and a = y/1 — + £. Notice 
that g{r s ,^) diverges logarithmically for £ — > 0, in agreement with our previous findings. 

Fig. 3 focuses on the paramagnetic degenerate regime: it presents a comparison between 
the numerical evaluation of from Eq. (2) (curve labelled C ) and its analytical approx- 
imation Eq. (12) (curve labelled B). The curve labelled as D is obtained neglecting the 
T 2 lnT term in Eq. (12): by comparing the D curve with the B curve, the importance of 
the T 2 InT correction becomes evident: without it the spin drag is strongly underestimated. 
The A curve represents the result obtained using in Eq. (2) the zero-temperature limit for 
ImQ 0a (q,ui). As expected, Fig. 3 shows that in this range of temperatures, the analytical 
approximation (B curve) is very close to such result. 



IV. THE NON-DEGENERATE LIMIT 

In the non-degenerate limit k B T » Ep, the dimensionless chemical potential fi a of Eq. 
(5) is given by the classical formula p, a = T a \n(h 2 27m a /m*k B T). Since —fL a /T a — > oo we 
can replace the function F(z) (Eq. (7)) by its limiting form F(z) — > ^fix /2z for z — > oo. 

This leads to the result 

T _ , x T>7> 1 E Fa / m*7r UJ — m *" 2 ,, fi V 

ImQ 0a (q,uj) w / e (16) 

a* fcsi V 2.KBJ- q 

t 2 2 
— n q 

identical to the 3D case. The Gaussian factor e 8m * fc s T on the right hand side of this equation 
tends to 1 in the classical limit (i.e., for h — > 0) for any finite wave vector. However, we 
cannot set H — 0, because the integral extends to arbitrarily large wave vectors, and the 
classical approximation fails at large enough values of q. The Gaussian factor assures the 
convergence of the integral in Eq. (2) by cutting off the integral at values of q such that 
h 2 q 2 /2m* ~ k B T. 
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Using the classical expression for the dielectric constant e(q,cu) ~ 1 + ko/q, with kr> = 
2ne 2 n/ekBT, after simple manipulations, we finally obtain 



t>>t f h 3/2 Ry* f°° y^e x 



,2 



where A = (2 3 / 2 /r 2 )(Ry* / K B Tf/ 2 and ify* = m*e 4 /2e 2 h 2 is the effective Rydberg. 

Remarkably, the last result (corresponding to k B T ^> Ry* or r s ^> 1) is independent of 
the electronic density. It differs from the analogous result in three dimensions, which scales 
as l/T 3 / 2 /n(T) and has a weak (logarithmic) dependence on the electronic density [6]. 

The limiting form of the transresistivity Eq. (18) is shown in Fig. 4 (dashed-dot line) along 
with the results obtained from Eq. (2) for the same choice of carrier densities as in Fig. 1. We 
stress that this curve delimits the region of the (n, T) plane occupied by the family of curves 
Pl I (n, T) and is approached at lower and lower temperatures for decreasing densities (see 
Fig. 4). The inset of Fig. 4 shows the comparison between pjj_ calculated for n = 10 l0 cm~ 2 
(solid line) and its limiting behaviors Eq. (17) (dashed line) and Eq. (18) (dashed-dot line). 
It is interesting to notice that, because of the low Fermi temperatures and at variance with 
the 3D case, in 2D the non-degenerate behavior of the spin-transresistivity could be observed 
experimentally, since it corresponds, for reasonable densities, to temperatures as low as 100K 
(see inset Fig. 4). We also notice that the non-degenerate limit Eq. (17) is in reality a good 
approximation even for temperatures as low as few Tp. 



V. DISCUSSION AND SUMMARY 

A few comments are now in order: as Eq. (12) shows, in the degenerate regime, the 
2D spin transresistivity behaves as T 2 (A + BlnT), in contrast to the 3D case in which 
PH ~ T 2 (see Ref. [4]). The additional logarithmic correction is due to the fast variation 
of ImQ 0rr (q,ui) for q fa 2kp a , so that the high q regime dominates the low-temperature 
behavior of p-^. This correction is also present, in principle, in the more familiar Coulomb 
drag transresistivity between two separate electron layers. However, in that case, due to 
the presence of the insulating barrier of thickness d, the Fourier transform of the interlayer 
Coulomb interaction is given by 27re 2 e~ qd / eq and decays exponentially for large q . Thus, in 
CD the influence of the logarithmic term on the behavior of p-|-j_ in the degenerate regime is 
negligible. 

The exponential decay of the interlayer interaction is responsible for another main differ- 
ence between the Coulomb drag and the spin Coulomb drag. As can be seen in Ref. [11], for 
T ~ Tp the transresistivity is dominated by the plasmon modes (especially the acoustic one), 
which enhance the contribution of small wave vectors. In the 2D-SCD, however, the acoustic 
mode, in which the two electron populations oscillate out of phase, is absent. Fig. 5 shows 
how the function I(u,q) = exp(— 2dq)ImQ ^(q, u)ImQo[(q, w)/4 \e(q,u)\ 2 smh. 2 (hu)/2kBT), 
behaves for small (upper panel) and intermediate (lower panel) values of q in the paramag- 
netic regime. For d — 0, q) corresponds to the integrand of Eq. (2); in the figure, how- 
ever, the interaction between the two spin populations is taken as v-^ = (2ire 2 /q) exp(-qd) 
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and the RPA dielectric function as e(q,u) = [1 — v(q)xo-\{q, w)] [1 — v(q)xoi(q, w)] — 
v uXo'\{ < 1j uj )Xoi{ ( 1j uj )- Different curves correspond to different values of d (0 < djkp < 9, as 
labelled). Let us first focus on the upper panel in which q = 0.1k F : for d — the main 
contribution to the integral comes from the single-pair excitation continuum. In this case the 
only possible collective mode (corresponding to the zeros of e(q,u>)) is the in-phase optical 
mode, whose contribution (the spike labelled as 0) is in any case negligible. As the expo- 
nential factor in the interaction is turned on (i.e. as the scattering rate between the two 
spin population is decreased), the strength of the integrand is transfered from the single-pair 
excitation continuum to the collective modes and the out-of-phase acoustic plasmon appears 
(see the two maxima labeled with 9, corresponding to d — 9kp)- 

The lower panel of Fig. 5 presents I(u),q) for q = 0.5/^: already at this intermediate 
value, the contribution to p^ becomes negligible when d > 0: in fact I(ou,q) scales as 

~ exp(— 2qd). 

We would like to emphasize a more general point: as Fig. 5 clearly shows, due to the 
absence of the exponential factor e~ qd in the Fourier transform of the Coulomb potential, 
the drag effect in the SCD is definitely larger than in the ordinary CD. In other words, in 
SCD, the two electron populations can transfer momentum to one another more effectively. 

The plasmon-enhancement is absent in 3D as well. In this case, however, this is due 
to a combination of three different effects: (i) the 3D plasmon frequency oj p {T) is finite at 
any temperature (u> p (T) > Vl67rna* 3 Ry*/h), (ii) due to the sinh 2 (/iu;/2£; B T) term in the 
denominator, the integrand in Eq. (2) is significantly large only for %uj ~ k B T, while it 
decreases exponentially for higher temperatures, and, (iii) the plasmon line-width, increases 
with temperature. Points (i) and (ii) imply that for small temperatures the integrand has 
already vanished when oo ~ uj p (T); on the other side, because of point (iii), for temperatures 
such that ksT ~ u p (T), the strength of the plasmon is negligible. 

Let us now examine some issues concerning the experimental observation of the SCD. 

First of all, it must be clear that the SCD is an intrinsic many-body effect, and, therefore, 
it is possible to design an experiment to measure it directly and independently of the ordinary 
diagonal part of the resistivity tensor. Such an experiment has been described in Refs. [4,6]. 

For different experimental setups, however, it might be important to know how the spin 
drag resistivity compares to the familiar Drude resistivity po- This information is provided 
in Fig. 6 for doped layers of InAs (m* = 0.026m e , e = 13.6) and GaAs (m* = 0.067m e , 
e = 12) at two different densities. Since it is customary to express the resistivity of doped 
semiconductor samples in terms of their mobility, in the upper panel we plot the quantity 
PH = —Xjnep^i as a function of temperature, pn has the dimensions of a mobility, and its 
value should be compared to that of the ordinary mobility po of the sample : if, at a given 
temperature, p^ < Pd the spin drag resistivity is larger than the ordinary Drude resistivity. 
It is evident from this figure that for appropriate but realistic parameters, p^ and pr> can be 
quite comparable. The significant numerical difference between the InAs and GaAs results 
is primarily due to the larger effective mass, and, to a lesser extent, to the smaller dielectric 
constant of GaAs. The impact of these two parameters on the spin drag resistivity is easily 
understood. Keeping in mind that the SCD is a consequence of Coulomb scattering between 
different spin populations, a smaller dielectric constant means stronger Coulomb scattering, 
and a larger effective mass implies a higher density of states available for scattering. Notice 
that the minimum value of p^ decreases with increasing density. 
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In the lower panel of Fig. 6, we plot the ratio P]J Pd (which is proportional to the sample 
mobility) versus temperature for a mobility po = 3 x 10 3 cm 2 /Vs. The same parameters 
and materials of the upper panel have been chosen. As the figure shows, for reasonable 
temperatures and realistic parameters, p^ can be a large fraction of pp. 

A peculiar signature of the two-dimensional SCD, which would be interesting to observe, 
is the non-degenerate density-independent behavior given by Eq.(18). As shown in the 
inset of Fig. 4, a sample would display such a behavior at T >> Tp, i.e. in the high 
temperature/low density [?] regime. Since, for densities of the order of n — 10 10 cm -2 , T F is 
as low as AK in GaAs and UK in InAs, this regime should be observable in these materials 
for temperature considerably lower than room temperature. 

Our calculations did not include the finite width of a realistic sample. The response of a 
quasi-2D sample would be somewhere in between our 3D and 2D results, depending on how 
well the idealized two dimensional approximation is fulfilled. In both the 2D and 3D cases 
however, we see that the ratio Ph/pd remains of the same order of magnitude for the same 
material and mobility and comparable density(see Fig. 6, lower panel, and Ref. [6], Fig. 2). 
Moreover even in the fully three dimensional case, the non-degenerate regime is only weakly 
(logarithmically) dependent on the carrier density. The above observations imply that a 
modest finite size correction would not modify significantly the relative importance of the 
predicted effect, or its non-degenerate behavior. 

In summary we have discussed the effect of dimensionality on the spin-transresistivity 
and shown that in two, as well as in three [5,6] dimensions, the spin Coulomb drag can 
be a sizeable effect. We have discussed the differences between Coulomb drag and spin 
Coulomb drag effects showing how, in the degenerate regime, the different form of the 
Fourier transform of the Coulomb interaction determines a quantitatively different behavior 
of the spin-transresistivity in comparison to the CD-transresistivity. We have also explained 
why in general the spin Coulomb drag is expected to be larger than the ordinary Coulomb 
drag. Finally we have demonstrated that, in the non-degenerate two-dimensional regime, 
the spin-transresistivity displays a universal density-independent behavior. Such a behavior 
could be in principle observed experimentally for realistic system parameters. 
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predictions on the spin Coulomb drag will also be qualitatively correct. 
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FIGURES 



FIG. 1. Spin transresistivity pn as a function of temperature (rescaled by Tp) for InAs param- 
eters (m* = 0.026m e , e = 13.6). Each curve corresponds to a different density as labeled. 

FIG. 2. Parameter space for the existence of a non-zero ImQoo-(q,uj) at T = 0. 
Upper panel: paramagnetic case {kp^/kp^ = 1). The three curves correspond to the set of points 
in which one or both of the arguments of the Heaviside functions in Eq. (8) are equal to zero: 
( ( //^f) 2 /2 + q/kp (solid curve), — (q/kp) 2 /2 + q/kp (dashed curve) and (q/kp) 2 /2 — q/kp (dotted 
curve). ImQo a (q,uj) is different from zero in the region between the solid and the dotted curve. 
The points labeled by A to D correspond respectively to qA a (^>) to (/Do-(w) for fmj/2Ep a = 0.2. 
Lower panel: spin-polarized case {kp\/kp\ = 1.5). The solid curve define the parameter space 
corresponding to ImQoi(q, w )' ^ ne dashed curve the one related to IrnQo^q, to). The points labeled 
as Bj and CI correspond to qsi and qc\ respectively, and delimit the interval relevant to the 
calculation of Eq. (14). 

FIG. 3. Degenerate regime: Comparison between p^ (curve labelled C), its analytical approxi- 
mation Eq. (12) (curve labelled B), the approximation obtained using the zero-temperature form of 
ImQo a (q,Lj) in Eq. (2) (curve labelled A), and the analytic one obtained by neglecting the T 2 InT 
term in Eq. (12) (curve labelled D). The curves are plotted vs temperature (rescaled by Tp) for 
n = W 10 cm~ 2 and InAs parameters. 

FIG. 4. Non-degenerate regime: Comparison between p^ (solid line) and its universal limiting 
behavior (dash-dotted line) Eq. (18) vs temperature for different carrier densities and InAs param- 
eters. Inset: Comparison between p-M (solid line) and its approximations (Eq. (17) (dashed line) 
and Eq. (18) (dashed-dotted line)). The curves are calculated for n = 10 10 cm~ 2 and plotted vs 
temperature (rescaled by the Fermi temperature). 

FIG. 5. Behavior of I(uj,q) = exp(—2dq)ImQo^(q,uj)ImQoi(q,u})/4\e(q,uj)\ 2 sinh 2 (h(jj/2kBT) 
vs the rescaled frequency tiuj/2Ep for fixed q and temperature but different values of d 
(0 < d/kp < 9). 

Upper panel: small-q 1 behavior (q = O.lfcp). Each curve corresponds to a different value of d/kp 
as labelled. 

Lower panel: intermediate-g behavior (q = 0.5kp). Curves are labelled with the corresponding 
value of djkp (as in the upper panel). 

FIG. 6. Upper panel: Spin drag mobility fj,^ = —1/p^ne vs temperature (rescaled by Tp) for 
n = W u cm~ 2 (dashed curve) and n = W 12 cm~ 2 (solid curve). Each couple of curves corresponds 
to a different material: InAs (m* = 0.026m e , e = 13.6) and GaAs (m* = 0.067m e , e = 12), as 
labelled. 

Lower panel: Ratio P\\/pd as a function of temperature for n = 10 n cm~ 2 (dashed curve) and 
n = W 12 cm~ 2 (solid curve) and sample mobility p = 3 x 10 3 cm 2 /Vs. Each couple of curves 
corresponds to a different material (InAs and GaAs , as labelled) 
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